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Abstract 

The surface detector (SD) array of the southern Pierre Auger Observa- 
tory will consist of a triangular grid of 1600 water Cherenkov tanks with 1.5 km 
spacing. For zenith angles 9 < 60° the primary energy can be estimated from 
the signal S(1000) at a distance of about 1000 m from the shower axis, solely on 
basis of SD data. A suitable lateral distribution function (LDF) S(r) is fitted to 
the signals recorded by the water tanks and used to quantify S(1000). Therefore, 
knowledge of the LDF is a fundamental requirement for determining the energy of 
the primary particle. The Engineering Array (EA), a prototype facility consisting 
of 32 tanks, has taken data continuously since late 2001. On the basis of selected 
experimental data and Monte Carlo simulations various preliminary LDFs are 
examined. 

1. Introduction 

High energy cosmic rays (CRs) are detected via the extensive air showers 
(EAS) they produce in the Earth atmosphere. Direction (9,<f)), energy (E) and 
mass of the primary CR are reconstructed from the secondary particles in the 
shower. The arrival times of shower particles at various detector locations give 
information on the arrival direction. The overall number of secondaries at ob- 
servation level scales roughly with primary energy, and the form of the shower 
depends to some extent on the primary mass. In the Auger experiment the lon- 
gitudinal shower development is measured by the Fluorescence Detectors (FD) 
while the lateral distribution at ground level is recorded by the SD, providing 
two independent measurements of the shower geometry and primary energy. The 
event reconstruction is hampered by the coarse sampling of the shower particles 
and by the statistical fluctuations of the shower development. High-developing 
showers are expected to have a flat lateral distribution, low-developing showers 
produce steeper lateral distributions. Fortunately, at about 1 km core distance 
the signal is virtually independent of primary mass and shower fluctuations, and 
is a good measure for the primary energy. Thus, the energy reconstruction re- 
quires, as a first and crucial step, to estimate S(1000) from few measured signals 
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at various distances from an a priori unknown core position. A second step is then 
to determine E from S(1000), which relies to a large extent on shower simulations 
and is therefore model dependent [1]. We do not discuss the energy calibration 
of the SD here. Since most of the Auger events have rather few SD stations hit, 
the reconstruction of the shower is not trivial. The functional form of an LDF, 
S(r), and its parameters, varies with 9, energy and mass, and its determination 
requires a good estimate of the core position, which in turn requires a reliable re- 
construction of the shower direction, which relies on a precise time measurement 
and stable trigger performance. The estimation of S(1000) (and thus E) can be 
greatly improved if the shape of the LDF is known. Here we present various 
approaches to determine the LDF from experimental data and MC simulations. 

2. Probing various LDFs 

In contrast to S(1000) the shape of the lateral distribution does not change 
much with energy. Therefore, it makes sense to decouple the normalisation con- 
stant from the shape parameters of an LDF and to combine showers of different 
energies. Simulations of EAS with AIRES/QGSJET01 in the range E = 1-100 
EeV and for 9 = 0-60° have been performed for the Auger experiment and their 
output was processed through a response simulation of detectors at core dis- 
tances 200-2500 m [2]. The LDF, in units of vertical equivalent muons (VEM), 
was parametrised with an empirical function of the form S(r) = E 0,95 • \q a + Bx + Cx 
with x = lg(r/1000 m) and the parameters A, B and C were determined as func- 
tion of 9 (see fig. El). Independent from simulations, the LDF was also deduced 
from experimental data. The EA was operating in very stable conditions during 
the period May to November 2002. Therefore data from this period have been used 
for the following analysis. High-quality events have been selected, which had a 
successful directional reconstruction with 9 < 60°, signals above 3 VEM in at least 
6 stations and a core position inside the EA. High-multiplicity events are very rare: 
only 3 x 10~ 3 of the events have 6 or more stations above threshold. Each event 
was examined and events or stations with obvious problems were removed from 
the sample. A few well defined events are better to determine the LDF than many 
events of lower quality. The following LDFs have been investigated: (i) a simple 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

sec 6 sec 9 sec 9 

Fig. 1. Shape parameters of LDF as function of sec 9 as predicted by MC simulations [2]. 
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Fig. 2. a) Measured LDF for two zenith angles (for power law fit). Full symbols: stations 
with signal, open symbols: silent stations. The data for 1 < scc6> < 1.25 have been shifted 
upwards for clarity. b) Residuals (S — S'th)/c r th as function of r for two values of 8 (for 
power law fit). Open symbols are silent stations. c) Fitted value of v as function of sec 9 
for individual events (squares) and averaged (stars). Solid line: fit to stars. d) Different 
LDFs for three zenith angles. For sec 8 = 1 a few error bars are plotted. 



power law S(r) = S'(IOOO) • (r/1000 m) v with a 9 dependent index v = a + bsec 8, 
(ii) an NKG-type function: S(r) = const, (r /r s )~ p ~ - 2 ■ (1 +r/r s )~ p = 5(1000) ■ 
(r/1000 m)^- a2 -((r + r s )/(1000 m + r,))-^ with (3 = a + bsecO and r s = 700 m*, 
and (iii) the MC inspired LDF S(r) = \0 A+Bx+Cx . These forms were fitted to 
individual events using a maximum likelihood fit of core location and LDF at the 
same time. Silent (i.e. alive but no signal above threshold) and saturated stations 
are properly included in the fit. The error cr(S)/S = yl).08 2 + 0.6/ S of a signal S 
(in VEM) is taken from an analysis [4] of data from a closely positioned detector 
pair. For the power law and the NKG-type LDFs two analyses were performed. 
First, in a two-parameter fit, the slope parameters v and (3, resp., have been var- 
ied together with the scale factor S(1000). Then a parameterisation of v and j3 
as function of 9 was determined, which was then used in a second analysis with 
only fitting S(1000). Fig. Et) shows the measured LDF (divided by S(1000)) for 



* Since (3 and r s are strongly correlated, we have fixed r s = 700 m and left (3 to vary. 
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Table 1. Moments (mean, a) of residual distribution of exp. data with various LDFs. Only 
events with > 6 stations were used in the present analysis. 



two zenith angles when the power law assumption is used for the core finding. To 
quantify the quality of the fit, residuals (S — S t h) / &th as a function of r are formed 
(see fig. Et>))- For a good LDF the residuals should scatter for all r symmetrically 
around with a variance of 1. Means and standard deviations of the distribu- 
tion of the residuals are used to compare different LDFs. For the power law fit 
the fitted values on v are shown as function of 9 in fig. Et)- It is evident that 
higher-multiplicity events have a smaller scatter. A line is fitted to the averages 
(stars) of the full symbols that yields v = 5.1(±0.4) - 1.4(±0.2) sec9. With the 
9 dependence of v fixed, the only fit parameter left for the LDF fit is S(1000), 
leading to more stable fit results for low multiplicity events. In the same way for 
the NKG function the variation of f3 was found: (3 = 3.3(±0.2) - 0.9(±0.2) sec6». 
Fig. I2~t0 shows that the three chosen LDFs, with suitably adapted parameters, 
agree well within the experimental errors. For completeness also the Haverah 
Park LDF [3] is shown, that predicts smaller densities at large core distances. 

3. Results 

The moments of the residual distributions from experimental data includ- 
ing silent stations obtained with different LDFs are listed in table El There is no 
major bias apparent and all the distributions have about the same widths, which 
indicates that all three LDFs describe well the data presently available. Previous 
experiments have shown that a pure power law cannot describe the shower signals 
at large core distances [3]. The results presented here are still preliminary as the 
statistics of events from the Auger EA, especially for energies > 10 19 eV, are small 
and since preliminary algorithms for directional and core reconstruction are used. 
In future each improvement on the statistics, angular and core resolution, and 
specifically hybrid events with their superior geometric reconstruction, will also 
improve the knowledge on the LDF and permit finer details to be analysed. 
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